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ABSTRACT 

Recent simulations of the formation of the first luminous objects in the universe predict isolated very 
massive stars to form in dark matter halos with virial temperatures large enough to allow significant 
amounts of molecular hydrogen to form. We construct a semi-analytic model based on the Press- 
Schechter formalism and calibrate the minimum halos mass that may form a primordial star with the 
results from extensive adaptive mesh refinement simulations. The model also includes star formation 
in objects with virial temperatures in excess of ten thousand Kelvin. The free parameters are tuned 
to match the optical depth measurements by the WMAP satellite. The models explicitly includes the 
negative feedback of the destruction of molecular hydrogen by a soft UV background which is computed 
self-consistently. We predict high redshift supernova rates as one of the most promising tools to test 
the current scenario of primordial star formation. The supernova rate from primordial stars peaks at 
redshifts 20. Using an analytic model for the luminosities of pair-instability supernovae we predict 
observable magnitudes and discuss possible observational strategies. Such supernovae would release 
enough metals corresponding to a uniform enrichment to a few hundred thousands of solar metalicity. 
If some of these stars produce gamma ray bursts our rates will be directly applicable to understanding 
the anticipated results from the SWIFT satellite. This study highlights the great potential for the 
James Webb space telescope in probing cosmic structure at redshifts greater than 20. 

Subject headings: cosmology: theory - early universe - supernovae: general 



1. INTRODUCTION 

The properties of pre-reionization luminous objects are 
integral to our comprehension of the process of reion- 
ization and their effect on subsequent structure forma- 
tion. Observations of distant {z = 6.28,6.4) quasars de- 
pict the relics of reioniz ation with their accompanying 
Gunn-Peterson troughs fB ecker et al.l 120021: iFan et aLl 
11002). Furthermore, Lyman alpha forest carbon abun- 
dances of ~10~^Zq and 1O~^'^Z0 observed at redshifts 3 
and 5, respectively, indicate that numerous early super- 
novae (SNe) enriched the intergalactic medium (IGM) 
(|Songaila fc CowidlToM iSongailal ImiH) . These early 
generations of stars are at least partly responsible for 
ionizing the Universe. A fraction of these stars lie in 
protogalaxies, but the other fraction of early stars are 
metal-free, form through molec ular hydrogen cooling and 
are very massive (M - IOOMq) l)Abel. Brvan. fc Normal 
12000. 2002). In ACDM cosmologies, these stars form at 
10 ^ z <, 50. The Wilkinson Microwave Anisotropy 
Probe (WMAP) data further constrains the epoch of 
reionization from the measurement of the optical depth 
due to electron scattering, t^s = 0.17 ± 0.04 (68%), which 
corresponds to a reionization redshif t of 17 ib 5 when as - 
suming instantaneous reionization l|Kogut et all 1200^ . 
However in this paper, we shall show that this epoch is 
gradual and its effect on primordial star formation. 

In hierarchical models of structure formation, small 
objects merge to form more massive structures. 
Eventually, a fraction of halos are massive 
5 X lO^ Mm at z '--^ 2 0) en ough to host cool- 
ing gas llCouchman fc Ree£lll986t iTegmark et alJli997t 
lAbel et a,].in998HFul]er fc Cnuchmanll200Hl . These halos 
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do not cool through atomic line cooling since T^ir < 10"*; 
however, primordial gas contains a trace of H2. Free 
electrons which allow molecular hydrogen to form that 
acts as an effective coolant at several hundred degrees 
through rotational and vibrational transitions. H2 is eas- 
ily photo-dissociated by p hotons in the Lyman- Werner 
(LW) ban d, 11.26-13.6eV jFieid. Somerville. fc DressleJ 
[1966; Stec her fc Williamj lT967jl . thus H9 can be de- 
stroyed by distant sources in a neutral Universe. Pri- 
mordial stars produce copious amounts of UV photons 
in the LW band, destroying the most effective cool- 
ing process at z ~ 20. This negative feedback from 
the UV background significantly inhibits the primordial 
star formation rate in the early universe by requiring a 
larger potential well for gas to condense. Many groups 
have explored the effects of a UV background on cool- 
ing and collapsing gas. Firstly, iDekel fc Reed l)1987fl 
discovered that H2 can be diss ociated from large dis- 
tances. Then more quantatively, 'Haima n. Rees. fc Loehl 
(1997) questioned how the collapse homogeneous, spher- 
ical cl ouds are affected by a UV ba ckground. In more de- 
tail, (Ifaima^^^e^^^^ei l)2000D determined whether a 
gas cloud collapsed by comparing the cooling time with 
the current lifetime of the cloud, which was calculated 
in the presence of solving the spherically symmetric ra- 
diative transfer equation along with time-dependent H2 
cooling functions. However, these studies only considered 
spherically symmetric cases while realistically these halos 
are overdensities within filaments. To combat this prob- 
lem, Machacek. Brvan. fc Abel (2001) employed a three- 
dimensional Eulerian adaptive mesh refinement (AMR) 
simulation to determine quantitative effects of a UV 
background on gas condensation. 
We consider lower mass stars to form in more massive 
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halos that may fragment via atomic line cooling as well as 
the metal lines from the heavy elements expelled by the 
earlier generation of Pop III star s. We use the prescrip- 
tion outlined in Haiman & Locb l)1997j) (hereafter HL97) 
to model the abundances and luminosities of these stars. 

With infrared space observatories, such as the 
Spitzer Space Tele scope, Primordial Explorer (PRIME; 
iZheng et al.ll200j) . and James Webb Space Telescope 
(JWST), sufficient sensitivity is available to detect the 
SNe from primordial stars. Although SNe remnants 
are bright for short periods of time, they may be the 
best chance to directly observe primordial stars due to 
their large intrinsic luminosities. If these events are 
recorded, many properties, such as mass, luminosity, 
metallicity, and redshift, of the progenitors can be cal- 
culated using metal-free or ultra metal-poor SNe mod- 
els, w hich can validate or falsify sim ulations of the first 
stars ijAbel. Brvan. fc Normanll2002(l . To evolve our cal- 
culation through redshift, we need to retain informa- 
tion about the UV background and number densities 
of primordial stars, which can extend to determining 
such quantities as volume-averaged metallicity and ion- 
ized fraction of the Universe. Our model is constrained 
with (a) the WMAP optical depth measurement, (b) pri- 
mordial star formation and suppression in high resolution 
hydrodynamical simulations, and (c) local observations 
of dwarf galaxies that constrain high redshift protogalax- 
ies properties. 

We organize the paper as follows. In §2, we report the 
semi-empirical method behind our calculations, which 
incorporates effects from negative feedback and an ion- 
ized fraction and is constrained from the WMAP result 
and local dwarf galaxy observations. In §3, we present 
the number density of primordial SNe in the sky and 
the evolution of SNe rates, metallicity, and related SNe 
quantities. We also explore the feasibility of observing 
these SNe by calculating their magnitudes and compar- 
ing them with the sensitivities of infrared space observa- 
tories. Then in §4, we suggest possible future observa- 
tions and numerical simulations that could further con- 
strain our model. Finally, §5 discusses and summarizes 
our results and the implications of the first observations 
of SNe from primordial stars. 

2. THE METHOD 

We use several theories and results from simulations 
of structure formation and metal-free stars. We use a 
ACDM cosmology with Q/, = 0.70, f^cDM = 0.26, fib = 
0.04, h — 0.7, (78 — 0.9, and n = 1. Qa, JIcdm, and 
fib are the fractions of mass-energy contained in vacuum 
energy, cold dark matter, and baryons, respectively, h is 
the Hubble parameter in units of 100 km s~^ Mpc~^. n = 
1 indicates that we used a scale- free power spectrum, and 
(Tg is the variance of random mass fluctuations in a sphere 
of radius 8h~^ Mpc. W e use a CDM power spectrum 
defined in iBardeen et a l. (1986), which depends on as 
and h. 

The remaining parameters are the primordial stellar 
mass, and the factors iesc and f^, which dictates the pro- 
duction and escape of ionizing photons, where fesc is the 
photon escape fraction, and f^, is the star formation effi- 
ciency. 

To evolve the ionization and star formation behaviors 
of the Universe, we must determine the density of dark 



matter halos that host early stars. In Tyir < 10^ K 
halos (henceforth "minihalos" ) , II2 cooling is the pri- 
mary mechanism that provides means of condensation 
into cold, dense objects. However in T^r > 10** K, which 
corresponds to My„ > 10*Mq[(1 + z)/10]-^^^, halos, 
hydrogen atomic line cooling allows the baryons to frag- 
ment and cool into stars. 

2.1. Minihalo Star Formation 

The first quantity we need to begin our calculation 
is the minimum halo mass that forms a cold, dense gas 
core due to H2 cooling. Radiation in the LW band photo- 
dissociates H2, which inhibits star formation in miniha- 
los. This negative feedback from a UV background does 
not necessarily prohibit the formation of primordial stars 
in minihalos, but only increases the critical halo mass in 
which condensation occurs, and delays the star forma- 
tion. From their simulations of pre-galactic structure 
formation, Machacck. Brvan. & Abel (2001) determined 
the minimum mass of halo that hosts a massive primor- 
dial star is 




. (1) 

where Mmin is the minimum halo mass that contains a 
cold, dense gas core, fed is the fraction of gas that is cold 
and dense, and F^w is the flux within the LW band in 
units of 10~^^ erg~^ s"^ cm^^ Hz^^. For our calculation, 
we consider fed = 0.02, which is a conservative estimate 
in which we have an adequate source of star forming gas. 
With our chosen fed and no UV background, a 1.74 x 10^ 
M(7) halo will form a cold, dense core, which will continue 
to form a primordial star. In a typical UV background of 
J = 10~^^ erg~^ s~^ cm~^ Hz~^ sr~^, the minimum mass 
is 4.16 X IO^Mq. Additionally, minihalos can only form 
a star within neutral region s of the Universe s ince they 
are easily photoioni zed (,Haiman. Abel. & Mad au 20Q2 
lOh fc HaimanI 12?)?)^ and H I is a necessary ingredient 
for producing H2. 

Nu merical simulations fe.g. |Al,)gl, Br van. fc NormarJ 
120021 iBromm. Coppi. fc LarsonI 12002 ^ illustrated that 
fragmentation within the inner molecular cloud does not 
occur and a single massive (M ~ IQQMq) star forms in 
the central regions. These stars produce hard spectra 
and tremendous amounts of ionizing photons. For ex- 
ample, the ionizing photon to stellar baryon ratio n-y 
91300, 56 700, and 5173 for H, He, Hc+ in a 2OOM0 star 
l|Schaereii r20021. We calculate the ionizing photon flux 
by considering 

where Q is the time-averaged photon flux, Tiife is 
the stellar lifetime, and Pmini is the comoving num- 
ber density of minihalos that we determine from an 
elliposdial variant of Press-Schechter (P S) formalism 
UPress fc SchechteHlT97ilsi;^h fc Tormen 20021. 

The primordial initial mass function (IMF) is un- 
known, therefore, we assume a fixed primordial stellar 
mass for each calculation. We run the model for primor- 
dial stellar masses, Mjs, of 100, 200, and 500 AIq and 
use the time-averaged emissivit ies from metal -free stel- 
lar models with no mass loss l|Schaeredl200^ . In our 
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models, the minihalo i s quickly ionized and all photon s 
escape into the IGM l|Whalen. Abel, fc Noriiianl 120041) . 

Furthermore, we use a blackbody spectrum at 10^ K to 
approximate the spectrum of the primordial star to con- 
sider the flux in the LW band since surface temperatures 
are virtually independent of mass at M >, SOM©. 

2.2. Star Formation in Protogalaxies 

In Tyir > lO"' K halos, baryons can cool efficiently 
through atomic line cooling, thus fragmenting and form- 
ing stars. As described in HL97, we parameterize the 
properties of these stars by a couple of factors. 

2.2.1. Star formation efficiency 

High redshift galaxies appear to have similar properties 
as local dwarf galaxies. We can constrain the star forma- 
tion efficiency f* by letting local observations gu ide us. 
In these galaxies, f^ range from 0.02 to 0.08 (Tavloi ^t al.l 
[1999: Walter et al. 2001). Using the orthodox Schmidt 
star formation law,iGnedin (2000) estimated f^, to be 0.04 
if it was constant over the first 3 Gyr but can also be as 
low as 0.022. It should be noted that f^, can be higher if 
star formation ceased after a shorter initial burst. 

Analyses of metallicities in local dwarf galaxies reveal 
their prior star formation. Both Type la and Type II 
SNe contribute iron to the IGM, but Type II SNe pro- 
vide most of the a-process elements (e.g. C, N, O, Mg) 
to the IGM. Type II SNe occur on timescales < 3 x 10'' 
yr wh ile Type la are delay ed by 3 x 10^ yr to a Hubble 
time (jMatteucci fc Recc hi "2001). Therefore, we expect 
an overabundance of a-process elements with respect to 
iron. [a/Fe] versus [Fe/H]^ plots help us inspect the evo- 
lution of the ISM/IGM metallicity. If the prior star for- 
mation is inefhcient (spirals, irregulars), [a/Fe] is only 
shortly overabundant, which is characterized by a short 
plateau versus [Fe/H]. On the other hand, if the star 
formation is fast and occurs early in the lifetime of the 
galaxy, [a/Fe] remains in the plateau longer due to the 
quick pr oduction of meta ls within Type II SNe (sec Fig- 
ure 1 in lMatteucc'ill2002|) . Recentlv. IVenn et all (2003) 
discovered no apparent plateau in [a/Fe] versus [Fe/H] 
in dwarf spheroidal and irregular galaxies. They con- 
clude that star formation must have been on timescales 
longer than Type la SNe enrichment, which hints at a 
low and continuous star formation rates in these proto- 
galaxies when compared to recent star formation. 

These low star formation efficiencies are further sup- 
ported by the galaxies contained in the Sloan Digital Sky 
Survey (SDSS; York^al. 2000). Star formation cfhcien- 
cies within low mass galaxies (M < 3 x lO^^M©) decline 
as M2/3 (,Kauffmann et al.ll200.'^ . At high redshift and 
before reionization, most of the protogalaxies tend to be 
few X 10^ Alp,, which is comparable to many local dwarf 
galaxies (Mateo 1998). It should also be noted that even 
within individual Local Group galaxies no star formation 
history is alike, but it is worthwhile to adopt a global 
star formation efficiency and observe the consequences 
on reionization and primordial star formation. With the 
stated constraints, we take f* = 0.04 i n our m ain model 
in concordance with the Schmidt Law l|Gnedinil20C)dj) and 
stellar abundances (|VeTin et aLl(200^ . 

^ We use the conventional notation, [X/H] = log(X/H) - 
log(X0/H0) 



TABLE 1 

Ionizing photon production per stellar baryon 

AND luminosities FOR METAL-POOR IMFS 
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"7, He + 


log C 










[erg s-i M-i] 


10-7 


20102 


8504 


8 


36.20 


10-5 


15670 


5768 


0.2 


36.16 


0.0004 


13369 


3900 





36.28'' 



"See text for a discussion on values. 



2.2.2. Stellar luminosities 

We considered several IMFs for star formation within 
protogalaxies. Our main model uses a Salpeter IMF 
with a slope a = -2.35, metallicity Z = 10-^, 10-^, and 
0.0004, and (M,,,,^, M„p) = (1, 100)Mq. We take the con- 
tinuous starburst spectrum evolution from thi s particu- 
lar IMF that was calculated with Starburst99 ijSchaered 
120031 ILeitherer et al.i r i9 99^. To estimate the ionizing 
photon per stellar baryon ratio, n^, and luminosities of 
the IMF at a particular metallicity, we interpolate in 
logio space between the two adjacent IMFs. When Z < 
10-^, then we consider the Z = 10"^ IMF. The proper- 
ties of these IMFs are listed in Table 1. Note that more 
metal-rich starbursts result in a softer spectrum and less 
ionizing photons in which n^^ h decreases almost by a 
factor of 2. 

In the Z = 0.0004 IMF, we retain the luminosity 
from the Z = 10"^ and set n..^^ He+ = 0. We choose 
to do so because of the theoretical uncertainty of hot 
Wolf-Rav et (WR) stars and their pr esence (for a re- 
view, see lSchaerejl2000l) . In lSchaered l|2003fl . the lumi- 
nosity and spectral hardness increases due to the pres- 
ence of hot Wolf-Rayet (WR) stars at higher metallici- 
ties. However, Smith. Norris. & Crowther (2002) calcu- 
late hot WR spectra that are significantly softer. 

2.2.3. Ionizing photon escape fraction 

Radiation emitted by these stars have a probability 
fesc to escape from the protogalaxy and ionize the IGM. 
The protogalaxy ISM density and comp osition plays the 
biggest role in determining this factor. iHeckman et alJ 
(2001) showed that local and distant starburst galax- 
ies, including the gravitationally lensed galaxy MS 1512- 
cB58 (z = 2.7), have fesc ^ 0.06. This result agrees 
with previous analyses of the fe^r. in the Lyman contin- 
uum (Leitherer et al."1995'; 'Hurwitz. Jelinskv. & Dixon| 
[1997). How ever, Ly man break galaxies may have higher 
f esc > 0-2 ffSteidel.'mtini. fc Adelber gor 2001). In con- 
junction with f^ = 0.04 and Mfs = [100, 200, 5OO]M0, 
the choice of fesc = [0.050, 0.033, 0.028], respectively, in 
our calculation results in the same optical depth from 
WMAP, Tes = 0.17, thus we use thes e values in the 
main models. iWood fc Loebl l|2000|) and lRicotti fc Shull 
l)200nt) have argued that the escape fraction may be very 
small due to the much higher densities at high redshift. 
However, the shallow potential wells and their small size 
make them also susceptibl e to photo-evap oration and ef- 
fects of radiation pressure ()Haehn elt 1995). Other uncer- 
tainities such as dust content, metallicity, and whether 
the ISM density scales as (l+z)-^ blur our intuition about 
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the escape of ionizing radiation. Therefore, we allow f^sc 
to vary from 0.001-0.25 since it is unclear whether high- 
redshift, low-mass protogalaxies allow photons to escape 
due to self-photoevaporation or absorb the photons due 
to a higher proper gas density when compared to star- 
burst galaxies. 

2.2.4. Ionizing Photon Rates 

These factors are multiplicative in amount of radiation 
that is available from protogalaxies to ionize the IGM. 
The rate of photons emitted that can ionize species X is 



~dt 



/, , N /esc /it ^7,X dipproto dz 

= po{l + z) '- ^ — , (3) 



proto 



/i nip 



dz dt 



where fi is the mean molecular weight, po — 
r26(3HQ/87rG), rup is the mass of a proton, and 'ipproto 
is the mass fraction contained in protogalaxies that is 
calculated by PS formalism. We restrict the product of 
fesc and fi, to be in a range from 10~''-10~^ since the 
resulting reionization histories fall within the measured 
Tes- The reionization epoch greatly depends on the fac- 
tors fesc and In order to explore the consequences of 
different values of t^s and their resulting SNe rates, we 
vary the factors fesc and f^,. 

2.3. Clumping Factor 

Overdense regions experience an increased recombina- 
tion rate. Overdensities are characterized by the gas 
clumping factor C = and the recombina- 

tion rate is increased by this factor. To calculate this 
parameter in our model, we ut ilize the same cosmolog- 
ical Eulcrian AMR code enzo | Brvan fc NormanlllQQflD 
as m fMachace k. Brvan. fc Aber (2001') in a 256 simu- 
lation with a comoving box side of 500kpc, eight levels 
of refinement, and the previously specified cosmological 
parameters. A total (gas) mass resolution of 1013 (135) 
Mq results from the resolution in the top grid, which 
is smaller than the cosmological Jeans mass Thus, 
we account for all collapsed halos in our clumping fac- 
tor calculation. The simulation is purely adiabatic with 
no background radiation or atomic/molecular cooling. 
We show Csim in Figure ^ Cooling only affects local- 
ized regions of star formation and does not contribute 



10- 



o 




Fig. 1. — The dashed line is the gas clumping factor as calculated 
in our adiabatic hydrodynamical AMR simulation. Using equation 
151 . we estimate the clumping factor in the ionized region to be 
solid line. 



greatly in the boosting of C. However, ionizing radiation 
causes phot oevaporation of halos, which de creases C in 
the process llHaiman. Abel, fc Madal]l200l . We correct 
for photoevaporation by considering all minihalos with 
Mj < M < Mmin are photoevaporated in the ionized 
regions of the Universe, where 



M, 



10" 



0.15 



-1/2 



0.02 



-3/5 



1 



11 



3/2 



(4) 

is the cosmological Jeans mass. For simplicity, we as- 
sume the IGM has a clumping factor of unity although 
underdense regions in the IGM correspond to Gjgm ^ 
1. We concentrate on the gas clumping factor in the ion- 
izing regions since the goal is to calculate the increase in 
recombination rates. Consider an ionized region whose 
volume filling fraction Fh is increasing. The expansion 
will incorporate unaffected, clumpy material into the re- 
gion. The first term in © accounts for this effect. Any 
overdensities in the ionized region will be photoevapo- 
rated in approximately a sound crossing time of the halo. 



ph 



Rv 



10 km s- 



(5) 



where /3 is a normalization factor that accounts for 
the differences in halo densities at various redshifts and 
masses and R^^^ is the mean virial radius of ha lo with 
Mj < M < Mmi„ ijHaiman. Abel, fc Madaul200H) . Since 
P remains within ^^15% of unity for minihalos, we infer P 
= 1. The photoevaporation of overdensities is the second 
term in the evolution of the gas clumping factor, which 
is 

C — 1 

C = CsimFu 7 • (6) 

'^ph 

The effect from photoevaporation is illustrated in Fig- 
ure ^ Overdensities that are engulfed by ionized regions 
cannot sufficiently increase the gas clumping factor to 
overcome the photoevaporation that occurs, and C re- 
mains within the range 1-3 for the entire calculation in 
our main model. This equation is weakly sensitive to tph 
since varying tph by a factor of 6 alters C only by a factor 
of 2. In a later paper, we shall computationally address 
the effect of radiation on the clumping factor. 

2.4. The Evolution of UV Background and Halo 
Densities 

We initialize the following method at z = 75 with no 
UV background, evolve the cosmological radiative trans- 
fer equation |(SJ at the questioned redshift, and repeat the 
described procedure until cosmological hydrogen reion- 
ization occurs. 

Given a minimum mass of a star-forming halo, we can 
exploit PS formalism to calculate number densities of 
these halos. For minihalos, we calculate their number 
densities for halos with masses above yimin O and virial 
temperatures below 10^ K. Likewise, the protogalaxy 
mass fraction is calculated by considering all halos more 
massive than a corresponding T^^^ — 10* K. 

We choose a variant of PS formalism, which is an 
ellipsoidal collapse model that i s fit to numerical sim- 
ulations (Sheth fc Tormcn"199;?: 'Shcth. Mo. fc TormenI 
[2001: Sheth fc Tormcn 2002) . This model is concisely 
summarized in FMo fc White! l)2002|) . In minihalos, it 




Redshift 

Fig. 2. — (Clockwise from upper right) Primordial SNe rates (yr"-*^ deg~^) per unit redshift; Cumulative primordial SNe rate (yr~^ 
deg~^); Specific intensity (erg~^ s~^ cm"'^ Hz~^ sr~^) in the LW band; Comoving density (Mpc~^) of halos above the critical star 
formation mass in neutral regions; Critical halo mass (Mq) for primordial star formation. The solid, dotted, and dashed lines correspond 
to calculations run with a fixed primordial stellar mass of 100, 200, and 500 Mq, respectively. 



is reasonable to assume one star forms per halo since 
Efcind < EsNe, where Ebmd and EsNe are the binding 
energy of the host halo and kinetic energy of the pri- 
mordial SNe, respectively. The gas is totally disrupted 
in the halo and requires ^lOOMyr ^ Ih, where Ih is 
the Hubble time, to re-collapse (Abel, Bryan, & Norman 
2002, unpublished). On main sequence, primordial stars 
evacuate la rge Stromgrcn spheres on the order of kpcs 
ijWhalen. Abel, fc N orman 2004,). Therefore, 

dz dt dz ' 

where is the comoving density of primordial stars. 
With that in mind, we can calculate the volume-averaged 
emissivi tv (II 111 of the U niverse by using metal-free stellar 
models (ISchaereI^l2002|) and the prescription outlined in 
§2.2 for minihalos and protogalaxies, respectively. 

We evolve the spectrum from early stars to investigate 
how the UV background behaves, particularly in the LW 
band, with the cosmological radiative transfer equation, 

^§-^^.Hiyj^-3HJ-c.J+^e, (8) 

where H — H{z) — \QQhE{z) is the Hubble parameter, 
and E{z) = + Vlm{\ + z)-^. J = J{v, z) is specific 

intensity in units of erg^^ s^^ cm^^ Hz^^ sr~^. k is the 



continuum absorption coefficient per unit length, and e = 
e{v] is the proper volume-averaged emissivity ijPeeblesI 
Il993^ . The emissivity will be the current total luminosity 
per unit volume of early stars, which can be expressed as 

e-,mni{v) = Pmmi X L{v) X f {v) X ( 1 - — ^ ) , (9) 

V riH J 

eproto{v) = X L{v, f-,) X f{v), (10) 

e{v) = (H- zf [eminiil") + iprotoif)] , (H) 

where f is a factor that accounts for absorption from 
Lyman series transitions. L is the luminosity of the ob- 
ject and is defined as 

[ 47ri?2B^(r ^ lO^X) (mini) 

^^"^ ^ ] fn f L BAT=2ZmQK) . ^ , 

y JOkJkJ- B^{„=2.7kT/hp;T=230aOK) \P'^'^^'~>) 

where is a blackbody spectrum, and C is listed in 
Table 1. R is the radius of the primordial star ijSchaereil 
l2002tl. fn ^ is the fraction of O stars in the starburst 
ljScha,ereHl20 03'l. k and hp are Boltzmann's constant and 
Planck's constant, respectively. For the protogalaxies, T 
^ 23000K because the spectrum in the LW band will be 
dominated by OB stars, and we weight the luminosity 
by this blackbody spectrum. This effect and the increas- 
ing UV background effectively squelches minihalo star 
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Photon Energy (eV) 

Fig. 3. — The red dashed Une is an example unprocessed spec- 
trum of a continuous source of radiation. After being absorbed 
and re-emitted by Lyman-series transitions, it transforms into the 
"sawtooth" spectrum {black solid line). This particular spectrum 
is for Zon = 30 and Zq^s = 20. 

formation. Emissivity will be nearly zero above 13.6eV 
in the neutral Universe due to the hydrogen and helium 
absorption. 

Photons from primordial stars and protogalaxies be- 
tween 13.6eV and several keV ionize the surrounding 
neutral medium. Using the intrinsic >13.6eV ionizing 
photon rates from primordial stars and protogalax- 
ies © , we calculate a volume-averaged neutral fraction. 
The change of ionized hydrogen comoving density due to 
photo-dissociation and recombination is 

where tihii and hh = 0.76po/mp ai'c the ionized and to- 
tal hydrogen comoving density, C is the clumping factor, 
and krec — 2.6 x 10^^^ s~^ cm~'^ is the case B recombi- 
nation rate of hydrogen at T w 10"*. 

In the case for helium, fc^ec is approximately equal to 
the hydrogen case, but since the helium number den- 
sity is less than hydrogen, less recombinations occur. 
Naively, this will result in a higher He+ fraction than 
H+ due to the hardness of the primordial radiation. Re- 
alisticly when the ionizing photons reach the ionization 
front, they will ionize either helium or hydrogen since 
hydrogen still has a finite photo-ionization cross-section 
above 23.6eV. We consider the He+ regions to be equal 
to the H+ regions. Then we add N^. He to N^^ h to com- 
pensate for this effect. We perform a similar analysis on 
the double ionization of helium in which k^ec, He+ =1.5 
X 10~^^ s~^ cm~^ but we evaluate equation H12I) directly 
with N^^ He+ ■ 

The absorption coefficient in Equation ((SJ is ignored 
since we take into account the Lyman series line absorp- 
tion by the following procedure. Photons with 11.26eV 
< E < 13.6eV escape into the IGM, which will photo- 
dissociate H2. To calculate the flux within the LW band, 
we must consider t he processing of photons by th e Lyman 
series transitions l|Haiman. Rees. fc Loeblll997|) . Before 
reionization, these transitions absorb all photons at their 
respective energies. These absorbers can be visualized as 
optically thick "screens" in redshift space, for which the 
photon must have been emitted after the farther wall 




Redshift 

Fig. 4. — {Top graph) The thick solid line is the volume- 
averaged metallicity from protogalaxies and primordial stars in 
the 2OOM0 model. The thin dotted, solid, and dashed lines are 
volume-averaged metallicities from 175, 200, 250 A/0 primordial 
stars only. {Bottom graph) Iron abundance relative to solar. Leg- 
end is the same as the top graph. The large spread in iron from 
primordial stars arises from the steep dependence of iron ejecta 
versus mass in pair instability SNe. 



in redshift. This process creates a sawtooth spectrum 
with minimums at redshift screens, which is illustrated 
in Figure|31 For instance, an observer at 2: = 20 observes 
photons at 12.5eV; it must have been emitted by the Ly7 
line at 12.75eV at z — 20.4. The fraction of photons that 
escape from these walls during an integration step is 

l-[(l + z)/(l + Z,,,.een)]^-^+" , . 

■' 1- [(l + z)/(l + z + Az)]i-5+"' ^ ' 

where a — -1.8 is the slope of a power-law spectrum 
(F^ (X v"") and 

^screen — (1 ~i~ ^obs) 1- (1^) 

VLyi is the nearest, blueward Lyman transition to v. In- 
herently, Zscreen must bc between z and z -\- Az. 

3. RESULTS 

In this Section, we present the results of our calculation 
of the evolution of SNe rates, volume-averaged metallic- 
ity, and UV background. We also present the variance of 
optical depth to electron scattering and SNe rates with 
different protogalaxy star formation scenarios. In Figure 
12 we plot minimum halo mass, density of those halos, 
metallicity, SNe rate, and UV background in the LW 
band for the main models. 

3.L Metallicity 

Metal-free stars with masses between ~14O-26OM0 re- 
sult in a SNe that is visible and eject heavy elements 
into the IGM. When a star is within this range, the 
stellar core has sufficient entropy after helium burning 
to create positron/electron pairs. These pairs convert 
the gas energy into mass while not greatly contribut- 
ing to pressure. This creates a major instability, where 
the star contracts rapidly until oxygen and silicon implo- 
sive burning occurs. Then the star totally disrupts it- 
self by these nuclear explosions ijBarkat. Rakavv. fc Sack! 
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Fig. 5. — Evolution of filling factors of ionized hydrogen (top) 
and doubly ionized helium (bottom). The legend is the same as 
Figure 121 



ll967HBond. Arnett. fc Carrll984(l . Stars between ~100- 
14OM0 experience this instabihty. It is not violent 
enough to disrupt the star, but pulsations and mass 
loss occur until equilibriu m is reached and t he hydro- 
gren envelope is e.icctcd ( Bond. Arnett. fc Ca rr 198^ 
iHeerer fc WooslevI 120 02') . Then since it creates a large 
iron core, it probably forms a black hole, preventing any 
further ejecta. Additionally, metal-free stars above M >, 
260Mq may result in direct black hole formation, which 
prohi bits any significant heavy element ejecta and after- 
glow l)Frver. Wooslev. fc Hegeill200U) . 

For primordial stars in the pair-instability SNe mass 
range, we determine a metallicity usin g our SNe densities 
and the metal production of these SNe ijHeeer fc WooslevI 
12002,) . Except for iron, significant spread in metallicity 
is not apparent for different choices of primordial stel- 
lar mass. Iron ejecta in pair instability SNe range from 
O.OO3-57M0, so naturally we expect a large spread. For 
star formation in protogalaxics, wc use carbon and metal 
production values from Rcnzini & Voh (198l| for 1-8 
Mq and lMaedeJ <fl992l|) for 9-120 M©. Additionally for 
15-25Mf7), we use the iron nucleosynthesis yields from 
lR,auscher et al.l l|2002D . At low metallicities and above 
~40Mq, direct black hole formation occurs after stel- 
lar collapse, which inhibits any ejecta from escaping the 
remnant. Furthermore above ~25M0, a black hole forms 
throu gh "fallback " of ejecta onto the neutron star rem- 
nant l|Frveil IT999j) . Thus we set an upper mass cutoff 
of 25M0 for metal enrichment from protogalaxics. With 
the chosen Salpeter IMF, 2.9 x 10-^ 0.023 and 0.071 Mq 
of iron, carbon, and metals, respectively, are produced 
for every Mq of stellar material. It should be noted that 
most of the carbon is produced by intermediate mass 
stars, and its injection is delayed by a few hundred mil- 
lion years, which corresponds to the first metal enriched, 
intermediate mass stars exploding at z ~ 12. Further- 
more, we ignore the delay time of any stellar lifetimes in 
FigureElto give a rough estimate of the volume-averaged 
metallicity evolution. 

With the evolution of metallicity, we can estimate the 
period of metal-poor (-5 <, [Fe/H] <, -1) star forma- 
tion. For the most metal deficient stars, their metal- 
licities of [Fe/H] — (-4.0, -5.3) correspond to a forma- 




0.01 0.1 

Time after SNe [yr (l+z)"'] 

Fig. 6. — Using the decay of ^^Ni and Equation I20i . we calculate 
the luminosity of primordial SNe and its effective temperature. 
Heavy (top) lines with circles are luminosities and light (bottom) 
lines are temperature. The blue dashed line wi th filled squares is a 
Type la luminosity evolution for comparison JWooslev &: Weaverl 
119861) . The lines from top to bottom (dotted, solid, dashed) are for 
stellar masses 175, 200, and 250 Mq. 



tion epoch of z ~ (28 , 36) fNorr is. Rvan. fc Beerjl2001l: 
[Christlicb ct al.l l2002D . Additionally, primordial stars in 
the 200Mq model account for [Fe/H] -5.4 of the to- 
tal iron ejecta. If we set the primordial stellar mass to 
be 25OM0, these stars eject tremendous amounts of iron 
into the IGM and possibly enhance [Fe/H] to -4.1. Con- 
versely, primordial stars at the lower mass range of pair- 
instability SNe produce insignificant amounts of iron and 
will not pollute the IGM. The IGM iron contribution 
from primordial stars is likely lower than these quoted 
values since in a proper IMF not all stars will result in a 
pair instability SNe. 

Possibly these metallicities correspond to the VMS 
(Very Massive Star) compon ent of metallicity of m ctal- 
poor (-4 < Z/Zq < -1) stars (|Qian fc Wasserburel2 002). 
They assume that a VMS contributes 1O~^Z0 to the 
IGM. By inspection of the evolution of metal yields from 
protogalaxics and primordial stars, the VMS component 
may contain contributions from protogalaxics since pri- 
mordial stars are suppressed by the UV background be- 
fore producing a universal metallicity of 10~^Zq. Lastly, 
metallicities found in z ~ 5 Lyman alpha clouds are 
[C/H] ^ -3.7 a nd much less than our calculated values 
l)Songailall2001j) . The bulk of the metals might not mi- 
grate to the void in which Lya clouds exist, and the 
host galaxies would retain the majority of the metals. 
This behavior may explain the lack of metals observed 
in Lya clouds and the higher than expected metallicity 
in the VMS component. However, the clustering of the 
first objects, primordial IMF, and the mixing of metals 
in the IGM will have to be taken into account before we 
can predict metallicities of the Lyman alpha forest and 
metal deficient stars with high confidence. 

3.2. Optical Depth to Electron Scattering 

Optical depth due to electron scattering H15|) is a good 
observational test to determine the neutral fraction of 
the Universe before reionization. The WMAP satellite 
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TABLE 2 

Gaussian fits to M^in for the main models 

log(M) = a X EXP[-(z-Zo)2/(2(t2)] + C 



Mq 


A 


Z() 


cr 


C 


100 
250 
500 


1.67 ± 0.01 

1.71 ± 0.00 

1.72 ± 0.01 


10.0 ± 0.2 
10.5 ± 0.0 
11.3 ± 0.2 


21.9 ± 0.1 
23.3 ± 0.0 
23.5 ± 0.1 


5.19 ± 0.01 
5.16 ± 0.00 
5.15 ± 0.00 



measured res(z) — 0.17 ± 0.04 at 68% confidence. 

Tesiz) = 



neCTTHC 



(15) 



where is the proper electron density and axH is the 
Thomson cross-section. To accurately calculate He, we 
must consider all ionizations of primordial gas, which in- 
cludes H+, He+, and He^"*". Therefore, the proper elec- 
tron density is 

rie = (1 + z)^ {nHFH+ + nHeFHe+ + 2nHeFHe++) , 

(16) 

where Fh+ F^g+ , and F^g++ are the ionized volume 
fraction of H+, He+, and He++, respectively. The effect 
of more luminous primordial stars is evident in Figure |31 
as they ionize the Universe faster at high redshifts. To 
match the WMAP result, less ionizing photons are neces- 
sary from protogalaxies if the primordial IMF is skewed 
toward higher masses. Although these models have the 
same total Tgs, cosmological reionization occurs at z = 
13.7, 11.6, 11.2 for Mps = 100, 200, and SOOMq, re- 
spectively. However, these reionization redshifts are not 
consiste nt with observed Gunn- Peterson tro ughs in z ^ 6 
quasars ("Beck er et al.l2002tlFan et alJ2002 ) and the high 
IGM temp eratu res a. t z ~ 4 inferred from Lya clouds 
ijHui fc Ha iman 200^ . Perhaps portions of the Universe 
recombine after complete reionizatio n, which wil l match 
the most distant quasar observations l)Cenl2003b(l . If this 
were true, the first reionization epoch has to be faster 
and earlier to compensate for this partial recombination 
and to match the WMAP result. This would lower our 
SNe rates slightly since primordial star formation will be 
further suppressed by the quicker reionization. 

The ionizing history of the Universe is directly related 
to the number of ionizing photons that are produced and 
escape into the IGM. If we fix t^s — 0.17, we constrain 
f* and fesc to a power law 



fesc — B 



(17) 



where B = [0.00307, 0.00239, 0.00217] and a = [0.906, 
0.839, 0.806] for primordial stellar masses 100, 200, and 
5OOM0, respectively. The flattening of the power law 
with increasing Mps indicates the increasing ionizing 
contribution from primordial stars. 

3.3. Primordial Supernovae Rates 

The natural units in our computation is comoving den- 
sity, yet a more useful unit is observed SNe yr~^ deg"'^. 
First we assume these SNe are bright for 1 yr and then 
correct for time dilation. We consider the equation 



d^N _ dVc dp^ 
dtdz dz dt 



i^ + z)- 



(18) 
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Fig. 7.— {Left to right) Magnitudes for Mps = 175, 200, 250 
Mq. The dashed horizontal lines indicate limiting magnitudes of 
space infrared observatories. The heavy {top) and light {bottom) 
lines are SNe at z = 15 and 30, respectively. The magnitudes for the 
175 Mq and z = 15 case is not shown since it is too dim. The solid, 
dotted, dashed, and dash-dotted lines are for Spitzer wavebands 
centered at (3.56, 4.51, 5.69, 7.96) /^m. Note: The detection limits 
of Spitzer and PRIME Medium are not shown since they are too 
high at 24.5 and 25.6, respectively. 



where the (1+z) ^ converts the proper SNe rate into the 
observer time frame, and 



dz 



D 



H- 



(1 



zYD\ 



E{z) 



(19) 



is the comoving volume element, il (deg^ — 3.046 x 10^** 
sr) is the solid angle of sky that we want to sample. = 
c/Hq is the Hubble distance. Da = DmI{^ + z) is the an- 
gular diameter distance, and Dm = D^ Jq E^^{z')dz' is 
the comoving distance ijPeeb lcs 1993). The above equa- 
tions are only valid for a flat ACDM universe. 

Since we only allow primordial stars to form in neutral 
regions, SNe rates from these stars are highly dependent 
on the fesc and f*, but less sensitive to the stellar pri- 
mordial mass, M^s. In our main models, the SNe rate 
varies little with primordial stellar mass and is 0.34 yr~^ 
deg~^. Even if we vary in a range 0.01-0.1 and flx Tg^, 
primordial SNe rates do not vary more than 10% from 
the main models when constrained by WMAP. As fgsc 
increases, the SNe rates decrease due higher ionized vol- 
ume filling factor. The effect of a higher f^ squelches SNe 
rates by two processes, a higher ionized filling factor and 
higher UV background, which forces primordial stars to 
form in more massive halos. 

The primordial SNe rate peaks at z ~ 12-20, later 
epochs for larger primordial stellar masses, and falls 
sharply afterwards due to the ensuing cosmological reion- 
ization. Primordial star formation ceases after z ~ 12-16. 
The combination of an increasing UV background, reion- 
izing Universe, and disruption of minihalos from primor- 
dial SNe suppresses all primordial star formation. For 
each main model, we fit a Gaussian curve to the mini- 
mum mass of a minihalo that can host a primordial star, 
and the parameters arc listed in Table 2 and are valid for 
z > zo. 

It is fascinating that some rare primordial SNe occur at 
z >, 40. As an exercise, we estimate the earliest epoch of 
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minihalo star formation in the visible Universe with PS 
formalism and by considering i t takes ^9.33 Myr for a 
halo to form a protostellar core l)Abel. Brvan. fc NormanI 
1^02). With PS formalism, the "first" epoch equals 
where the minihalo density is the inverse of the comoving 
volume inside z = 1000 (10523 Gpc'^). Then we include 
an additional 9.33 Myr for the ensuing star formation. 
The halo masses of 1.74 x 10^ and 10^ M© correspond 
to formation times of z ~ 71 and 64, respectively. An- 
other interesting event to calculate is when the SNe rate 
equals one per sky per year, which occurs at z ^ 51 when 
considering the collapse and sta r formation timescales . 
This epoch is in agreement with iMiralda-Escudel l)2003l ^ 
estimate of z ~ 48. 

3.4. Magnitudes and Observability 

The magnitudes of primordial SNe are important as 
their occurrence rates to catch these events unfolding in 
the distant universe. We exploit the ^^Ni output from 
metal-free SNe models to calculate lum inosities L(t) from 
the t wo-step decay of ^^Ni to ^^Fe ijHeger fc WooslevI 
l2r)02|l . We consider the emission spectrum to be a black- 
body spectrum with a time-dependent temperature H2U|I 
calculated using free expansion arguments. 



T{t) = 



L{t)M 



-,1/4 



SiraEt^ 



(1 



(20) 



where M is the stellar mass and E is the kinetic en- 
ergy of the SNe. The temperature and luminosity of 
a 175, 200, and 250Mo primordial star SNe are depicted 
in FigureEl Kinetic ener gy is taken from the SNe models 
lIHeger fc Wooslevll200l . 

In the upper range (M > 2OOM0) of pair- 
instability SNe, thes e events produce 1-57Mq of ^^Ni 
(jllcger fc Wooslevl2002i') . which will produce tremendous 
luminosities comp ared to typical Type II SN e outputs of 
only 0.1-0.4 Mq (|Wooslev fc Weaved IT986|) . However, 
uncertainty in the primordial IMF places doubt in the 
frequency of primordial star formation with high ^^Ni 
yields. Using conventional Type II SNe parameters of 
L ~ 3 X 10^2 gj.g g-i^ T ^ 25000K (2 days < t < 7 
days) , and T = 7000K (7 days < t < 2 months) that are 
tuned by observed light curves. iMiralda-Escude fc Rees I 
ljl997l) determine apparent magnitudes that are 1-3 mag 
lower than our values, which is due to the greater nickel 
production of new pair- instability SNe models. Their 
model should apply to the pair-instability SNe with lit- 
tle ^^Ni ejecta. The temperatures are higher than our 
values because we allow for free expansion of the ejecta 
due to a lower external medium gas density. Radiation 
from the primordial star should expel most of surround- 
ing medium to create a low density, highly ionized region 
of approximately 100 pc in siz e for a 120Mp, and bOOMp, 
primordial star, respectively ijWhalen. Abel, fc NormanI 
|2004). 

In typical Type II light curves, radioactive decay does 
not significantl y contribute to the luminosity in the 
plateau stage |PodovI I1993D . In primordial SNe, the 
^^Ni decay may overwhelm the typical sources of energy 
within the expanding fireball and create a totally differ- 
ent light curve. We stress that our light curve is a very 
rough estimate at the processes occurring within a pri- 
mordial SNe. 




-50 ! 



J 40 



Fig. 8. — The solid line is the mean proper separation between 
minihalos. The dashed line is the redshift where minihalos first 
receive radiation from neighboring sources forming at the same 
epoch. 



In Figure 13 we compare maximum AB magnitudes 
of low- and high-mass pair-instability SNe at various 
redshifts and sensitivities of space infrared observato- 
ries. We consider the sensitivities of SIRTF at 3.5^m, 
PRIME medium and deep surveys, and JWST with ex- 
posure times of 100s and 10^ s. We also give the typical 
light curve for a Supernova Type la. Clearly our simple 
estimate leads to very high predicted luminosities and 
more detailed numerical models of these explosions are 
clearly desirable. 

3.5. Other factors 

Other factors could alter the feasibility of observ- 
ing primordial SNe. For instance, ~2% of primordial 
stars might die in coUapsar gamma-ray bursts (GRBs) 
|Hegcr ct al. 2002). Lamb (1999) also provides a ratio 
of GRBs to Type Ib/c SNe rates that is ~3 x 10"^ 
{fbeam / 10"^)""^. Howcver, this estimate is for a nor- 
mal Salpeter IMF, and as mentioned before, the pri- 
mordial star IMF could be skewed toward high masses, 
which would increase the probability of a GRB. Zero- and 
low-metallicity, massive stars outside the pair-instability 
mass range die as GRBs or jet-driven SNe, which are sim- 
ilar to GRBs but not as energetic and spectrally hard. If 
we consider a proportionality constant and beaming fac- 
tor, we estimate that the all-sky primordial GRB rate. 



Rgrb 



0.01 



Igbb 
0.02 



GRBs yr" 



(21) 



In M ;< 140Mq stars, a collapsar results when it 
forms a proto-neutron stars, cannot launch a super- 
nova shock, and directl y collapses to form a black 
hole in ~1 s (iWooslevI 119931 iMacFadven fc WooslevI 
11999(1 . In M > 26OM0 stars, it does not create a 
proto-neutr on star and direct ly collap ses i nto a massive 
black hole ijFrver. Wooslev. fc Hegeni200ll) . Afterwards, 
these black holes can accrete gas a nd produce X-rays 
that can further ionize the Universe fRicotti fc OstrikeiJ 
J|)04b). A few of these X-r ay sources could be detected 
in th e Chandra deep f i eld (IRicotti. Ostriker. fc Gnedml 
I2OOI [Alexander eFaJ] I2003D. In m ost GRB models 
(for an overview, see lMeszaro3l2002|) . the radiation is 
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beamed to small opening angles due to relativistic ef- 
fects, which would render a fraction of GRBs to be unob- 
servable in our perspective. Finally, gravitational lens- 
ing will significantly increase the apparent magnitudes 
for selected primordial SNe for very small survey areas. 
However, the overall magnitude distribution is slightly 
dimmed by gravitational lensing l|Marri fc FerraralllQQSt 
iMarri. Ferrara. fc Pozzettill200o(l . 

Using the proper minihalo density, we calculate the 
average light travel time between sources, which indi- 
cates when negative H2 feedback first affects star forma- 
tion. Star formation occurs ~9.33 Myr after halo forma- 
tion, so radiation escapes into the IGM at a time trad = 
tH(zform) + 9.33 Myr, where Zform is the halo formation 
redshift. These sources have a mean proper separation 

— 1/3 

of d = / (l-|-z); therefore the time where radiation 

influences other halos is trad + d/c. The minimum of this 
time, corresponding to z = 54, is where radiation reaches 
other halos for the first time. The mean proper separa- 
tion and epoch of first radiation effects are illustrated in 
Figure 

Deviations in the amp litudes of fluctuations, asi a nd a 
running spectral model ijLiddle fc Lvthlll992l ll993l also 
affect primodial SNe rates. We ran a set of models with 
(Js = 0.8, and predictably, the rates decreased by a factor 
of ^2.5 due to lesser powers at small mass scales. When 
we match r — 0.17, rates are 0.12-0.2 yr^^ deg^^ with 
the rates peaking earlier at z ~ 20, but now the red- 
shift peaks are nearly independent of redshift with only 
(5z « 1 separating the 100 and 500Mq models. If we 
keep the protogalaxy parameters from the main mod- 
els, rates do not change; however, the rates peak later 
than our main models by Sz « 4, and r lowers to 0.14. 
According to inflationary models, the spectral index of 
flucutations should be sl owly varying with scale. While 
analyzing WMAP data, IPeiris et al.l lj2003 '^ determined 
that the flucutation amplitude is significantly lower at 
small scales. With lesser powers at small scales, primor- 
dial SNe rates decrease such as the case of a lower erg. 
In principle, studying primordial SNe occurrences with 
respect to redshift could furnish direct constraints on the 
power spectrum at these small scales. 



4. CONSTRAINING THE FREE PARAMETERS 

We have demonstrated our model's dependence on the 
free parameters, /eso /*) and Mfg. We suggest several 
observational and numerical methods in order to con- 
strain our models. 



4.1. Observational 

To constrain early star formation history of dwarf ir- 
regulars with the [a/Fe] vs. [Fe/H] comparison, metal 
analyses of stars with [Fe/H] < -1.5 are needed to de- 
termine early star formation rates o f particular systems 
(See Figure 1.5 in lVenn et alJl200'3l) . Dwarf spheroidals 
exhibit greater variance in star formation histories from 
system to system, but a similar study will increase our 
knowledge of high redshift star formation in these small 
galaxies to further constrain fi, at high redshifts. To 
improve on these studies, the advances in multi-object 



spectrographs (e.g. VIMOS^, FLAMES^, and CMOS'') 
enable the gathering of many stellar spectra in one ex- 
posure. This will greatly increase the stellar population 
data of dwarf galaxies. 

Since some massive stars die as a long-duration GRB, 
these events can convey information from the death of 
both Pop II and HI stars. The propagation of the ini- 
tial burst and afterglow provide information about the 
total energy, gas density in the vicinity, and the Lorentz 
factor of the beam. The host galaxy ISM properties will 
help constrain the fesc i n high redshift gal axies. As in 
the case of GRB 030329 l|Stanek et al.ll2003() . the power- 
law spectra of the afterglow can be subtracted to obtain 
a residual that resembles a typical SN spectrum, which 
may be used to roughly determine the mass of the pro- 
genitor. Observing the afterglow is necessary to deter- 
mine its redshift. The prospect of observing prompt af- 
terglows will be accomplished easier with SwiftP, which 
can possibly detect GRB afterglows to z = 16 and 3 3 
in the K and M bands, respectively ijGou et al.ll2004j) . 
Furthermore, > 50% of GRBs occur earlier than z = 5, 
and 1 5% of those high redsh ift GRBs are detectable by 
Swift ^ romm fc Loebll2?)0l . The comparison to nearby 
SNe can provide crucial information of the high redshift 
ISM. Finally, it is a possiblity to explore the intervening 
absorption with the fast-pointing a nd multi-w avelength 
observations of Swif t (Vrecswiik. Mo Uer. fc Fv nbo 200^ 
lLoebl2003trBarkana fc Loebi2004,l . This IGM absorption 
would constrain the reionization history of the Universe 
better, which may change the primordial SNe rates, but 
more specifically the rate per unit redshift which would 
roughly conform to the filling factor evolution. 

The radiation from protogalaxies and primor- 
dial stars will not only ionize the Universe but 
also contribute to the near-infrared background 
(NIRB). Calculations have shown that primordial 
stars can pro vide a significant fraction of radiation 
to the NIRB llSantos. Bromm. fc Kamionkowskil I200I 
iSalvaterra fc Ferrarall2003j) : however, the paradigm of 
early reionization set by WMAP was not considered at 
the time, f^ of protogalaxies and densities of primordial 
stars can be further constrained if future studies of 
the NIRB consider the large emissivities of zero- and 
low-metallicity sources at z ^ 15. 

Other outlooks include detecting high redshift radio 
source s and searching for 21cm absorption and emis- 
sion (Hogan fc Rccs 1979; Scott & Rccs 1990:' iliev et alJ 
1200 3: Furlanetto. Sokasian. fc Hernauist 2004) . By look- 
ing for 21cm signatures, observations would be directly 
probing the neutral regions of the Universe since the 
Gunn-Pcterson trough saturates only at a neutral filling 
factor of ^10^^. Lastly, additional sea rches for hig h red- 
shift starbursts fe.g. .Elhs et alJl2001l: IPello et alJl2004|) 
in lensed fields will furnish an understanding of the char- 
acteristics of these objects and help tighten our models 
of early "normal" and primordial star formation. 



4.2. Numerical Simulations 

^ http:/ /www. eso.org/instruments/vimos/ 
^ http:/ /www. eso.org/instruments/flames/ 

http: / /www. gemini.edu/sciops/instrumcnts/gmos/gmosIndex. html 
^ http://swift.gsfc.nasa.gov/ 
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Metal production from primordial stars cannot account 
for the volume filling fact or of metals as seen in Lya 
clouds in current models ijNorman. O'Shea. fc PaschosI 
1^004). Therefore, ubiquitous star formation in pro- 
togalaxies most likely polluted the sparse regions of 
the Universe. Combining the primordial and proto- 
galaxy metal output, simulations with proper metal 
transport should agree with the abundances observed in 
Lya clouds. Such simulations that takes into account 
both of these m etal sources and mat ches the results of 
iSongaila Sz Cow ic (199f|) and iSongaila ( 2001 ) is needed 
to constrain star formation before reionization. 

Numerical simulations are also needed to investigate 
radiative transfer from a protogalaxy stellar popula- 
tion. This scenario will contain more complexities 
than a single primordial s tar in a spherical halo as in 
iWhalen. Abel, fc NormanI l)2004|) . However a protogalac- 
tic simulation with the Jeans length resolved, star forma- 
tion and feedback, and radiative transfer, it will be pos- 
sible to study the evolution of the ISM in a protogalaxy, 
which will produce insight and better constraints on the 
photon escape fract ion, .fp.sn in the Lyman c ontinuum 
at high redshifts. iRicotti. Gnedin. & Skill l|2n02albr 
thoroughly study radiative transfer around protogalax- 
ies; however, f^sc is still a parameter when it should 
be determined from the radiative transfer results in the 
simulations. Ideally, the analytica l ideas about ion- 
izing t he Universe in iMadal lll995D: [Haardt k Madaul 
l|199()D : iMadau. Haardt. fc ReeA l|1999D should be real- 
ized in such simulations but in the context of the cur- 
rent paradigm of a high rcdshift reionization as indi- 
cated by WMAP. Also we may hope by using numeri- 
cal radiative transfer techniques for line and continuum 
radiation in three dimen s ion to push the simulations of 
lAbel. Brvan. fc NormanI 1)2002(1 to follow the entire ac- 
cretion phase of the first stars. This would lead to 
stronger constraints on the masses of the very first stars. 

5. SUMMARY AND CONCLUSIONS 

The WMAP measurement of optical depth to elec- 
tron scattering places a constrain of early cosmologi- 
cal reionization of the universe, which we show to be 
mainly from star formation in protogalaxies. This gen- 
eral result is in agreeme nt with othe r studies of reioniza- 
tion after WMAP (e.g. "Cen '2003a': 'SomcrviUe & Livi^ 

Foo'^at]So^sia^™a^'ll'^^ ^'^^^ -^^^"^^^ ^ Ostrike: 



primordial stellar masses. These results are upper 
limits since the rate of visible primordial SNe de- 
pends on the IMF because only a fraction will lie 
in the pair instability SNe mass range. The other 
massive primordial stars might result in jet-driven 
SNe or long duration GRBs. 

• Stellar abundances and star formation in local 
dwarf galaxies aid in estimating protogalaxy char- 
acteristics. We choose f* = 0.04 and iesc — 0.050 
in our 100 Mq model. In protogalaxies, the star 
formation efficiency is slightly lower than local val- 
ues, but the photon escape fraction is within local 
observed fractions. 

• Primordial stars enrich the IGM to a maximum 
volume-averaged [Fe/H] ~ -4.1 if the IMF is skewed 
toward the pair-instability upper mass range. A 
proper IMF will lessen this volume-averaged metal- 
licity since only a fraction of stars will exist in this 
mass range. If complete metal mixing does not oc- 
cur, there will be metal-rich and metal-poor regions 
relative to the universal metallicity. 

• The entire error bar of the WMAP measurement 
of optical depth to electrons can be explained by 
a higher/lower f^, and iesc- Protogalaxies can ion- 
ize the IGM easily since low metallicity starburst 
models produce 20-80% more ionizing photons as 
previously used (Z — 0.001) IMFs. Massive primor- 
dial stars provide ^^10% of the necessary photons 
to achieve reionization. No exotic processes or ob- 
jects are necessary. 

• Only the upper mass range of pair instability SNe 
will be observable with JWST since the low mass 
counterparts do not produce enough ^^Ni to be very 
luminous. 

Although the IMF of primordial stars is unknown, sim- 
ulations have hinted that a fraction of metal-free stars 
may exist in the pair-instability mass range. When ob- 
servational rates, light curves, and spectra are obtained 
from future surveys, these data would provide very strin- 
gent constraints on the underlying CDM theory as well 
as our understanding of primordial star formation. 



The radiation from protogalaxies squelches out pri- 
mordial star formation, and primordial SNe 
deg~^ yr~^ are expected. SNe rates can vary from 
0.1 to >1.5 deg~^ yr~^ depending on the choice 
of primordial stellar mass and protogalaxy param- 
eters while still constrained by the WMAP result. 
The peak of SNe rate occurs later with increasing 
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